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Abstract 

A measurement of the cross section of top quark pair production in proton-proton collisions 
recorded with the ATLAS detector at the LHC at a centre-of-mass energy of 7 TeV is reported. The 
data sample used corresponds to an integrated luminosity of 2.05 fhr 1 . Events with an isolated elec- 
tron or muon and a t lepton decaying hadronically are used. In addition, a large missing transverse 
momentum and two or more energetic jets are required. At least one of the jets must be identified as 
originating from a b quark. The measured cross section, <x rJ = 186 ± 13 (stat.) ± 20 (syst.) ± 7(lumi.) pb, 
is in good agreement with the Standard Model prediction. 
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Abstract 

A measurement of the cross section of top quark pair production in proton-proton collisions recorded with the ATLAS detector 
at the LHC at a centre-of-mass energy of 7 TeV is reported. The data sample used corresponds to an integrated luminosity of 
2.05 fbT 1 . Events with an isolated electron or muon and a r lepton decaying hadronically are used. In addition, a large missing 
transverse momentum and two or more energetic jets are required. At least one of the jets must be identified as originating from a b 
quark. The measured cross section, cr t - t = 186 ± 13 (stat.) + 20 (syst.) + 7(lumi.) pb, is in good agreement with the Standard Model 
prediction. 
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1. Introduction 

Measuring the top quark pair (it) production cross section 
(a,-;) in different decay channels is of interest because it can 
open a window to physics beyond the Standard Model (SM). 
In the SM, the top quark decays with a branching ratio close to 
100% into a W boson and a b quark, and it pairs are identified 
by either the hadronic or leptonic decays of the W bosons and 
the presence of additional jets. ATLAS has previously used the 
single-lepton channel JJ], and the dilepton channels including 
only electrons and muons [2] to perform cross-section measure- 
ments. 

The large cross section for it production at the LHC provides 
an opportunity to measure o~ t - t using final states with an elec- 
tron or a muon and a r lepton with high precision. The cr t - t in 
this channel has been measured at the Tevatron with 25% pre- 
cision J3] and recently by the CMS Collaboration at the LHC 
with 18% precision [7]. A deviation from cr t - t measured in other 
channels would be an indication of non-Standard Model decays 
of the top quark, such as a decay to a charged Higgs (H + ) and 
a b quark with H + decaying to a t lepton and a t neutrino, 
or contributions from non-Standard Model processes (Q 0]. 
ATLAS has set upper limits on the branching ratio of top quark 
decays to an H + bosons decaying to a t lepton and a neutrino 

a. 

This analysis uses 2.05 fb 1 of data collected by ATLAS at 
the LHC from pp collisions at a centre-of-mass energy of 7 TeV 
between March and August 2011. After application of kine- 
matic selection criteria that require one top quark to decay via 
W — > tv (where I is either a muon or an electron) and identifica- 
tion of a jet as originating from a b quark (b-tag), the dominant 
background to the it — > I + t + X channels with the r lepton 
decaying hadronically is the it — > I + jets channel in which a 
jet is misidentified as a hadronic t lepton decay. Therefore, r 



lepton identification (r ID) is critical for separating signal and 
background. The r ID methodology employed in this analysis 
exploits a multivariate technique to build a discriminant |0]. 
A boosted decision tree (BDT) algorithm is used The 
number of r leptons in a sample is extracted by fitting the dis- 
tributions of BDT outputs to background and signal templates. 
The results are checked using an alternative method, referred to 
as the "matrix method", based on a cut on the BDT output. 

2. ATLAS Detector 

The ATLAS detector 111 ill at the LHC covers nearly the en- 
tire solid angle around the collision point. 1 It consists of an 
inner tracking detector surrounded by a thin superconducting 
solenoid, electromagnetic (EM) and hadronic calorimeters, and 
an external muon spectrometer incorporating three large super- 
conducting toroid magnet assemblies. The inner tracking de- 
tector provides tracking information in a pseudorapidity range 
|?7l < 2.5. The liquid-argon (LAr) EM sampling calorimeters 
cover a range of I77I < 3.2 with fine granularity. An iron- 
scintillator tile calorimeter provides hadronic energy measure- 
ments in the central rapidity range (\t]\ < 1.7). The endcap 
and forward regions are instrumented with LAr calorimeters for 
both EM and hadronic energy measurements covering I77I < 4.9. 
The muon spectrometer provides precise tracking information 
in a range of I77I < 2.7. 



1 Atlas uses a right-handed coordinate system with its origin at the nominal 
interaction point in the centre of the detector and the z-axis along the beam 
pipe. The x-axis points to the centre of the LHC ring, and the y-axis points 
upwards. The azimuthal angle (p is measured around the beam axis and the 
polar angle 8 is the angle from the beam axis. The pseudorapidity is defined 
as rj = - ln[tan(#/2)]. The distance AR in r\ - <p space is defined as AR = 
V(A0) 2 + (A7,)2. 
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ATLAS uses a three-level trigger system to select events. The 
level- 1 trigger is implemented in hardware using a subset of de- 
tector information to reduce the event rate below 75 kHz. This 
is followed by two software based-trigger levels, level-2 and 
the event filter, which together reduce the event rate to about 
300 Hz recorded for analysis. 

3. Simulated Event Samples 

Monte Carlo (MC) simulation samples are used to optimise 
selection procedures, to calculate the signal acceptance and to 
evaluate contributions from some background processes. 

For the it and single top-quark final states, the next-to- 
leading-order (NLO) generator MC@NLO U2i is used with a 
top-quark mass of 172.5 GeV and with the NLO parton distribu- 
tion function (PDF) set CTEQ6.6 1 1 311 - The "diagram removal 
scheme" is used to remove overlaps between the single top- 
quark and the it final states. The it cross section is normalised 
to the prediction of HATHOR (164+ Jg pb) 0], which employs 
an approximate next-to-next-to-leading-order (NNLO) pertur- 
bative QCD calculation. 

For the background channels, MC samples of W/Z, single 
top-quark events and diboson WW, WZ, and ZZ events (all in 
association with jets) are used. W+jets events and Z/y*+jets 
events (with dilepton invariant mass ni(+(- > 40 GeV) are gen- 
erated by the ALPGEN generator [15] with up to five outgoing 
partons from the hard scattering process, in addition to the vec- 
tor bosons. 2 The MLM matching scheme of the ALPGEN 
generator is used to remove overlaps between matrix-element 
and parton-shower products. Parton evolution and hadronisa- 
tion is handled by HERWIG 1 16], as is the generation of dibo- 
son events. The leading-order PDF set CTEQ6L is used for all 
backgrounds described above. 

All samples that use HERWIG for parton shower evolution 
and hadronisation rely on JIMMY [ 17] for the underlying event 
model. The r-lepton decays are handled by TAUOLA [18]. 
The effect of multiple pp interactions per bunch crossing ("pile- 
up") is modelled by overlaying simulated minimum bias events 
over the original hard-scattering event [19]. MC events are then 
reweighted so that the distribution of interactions per crossing 
in the MC simulation matches that observed in data. The aver- 
age number of pile-up events in the sample is 6.3. After event 
generation, all samples are processed with the GEANT4 lEoll 
simulation of the ATLAS detector, the trigger simulation and 
are then subject to the same reconstruction algorithms as the 
data li2lll. 



4. Data and Event Selection 

The event selection uses the same object definition as in the it 
cross-section measurement in the dilepton channel Q with the 
exception of a r candidate instead of a second electron or muon 
candidate and some minor adjustments. The electrons must be 



isolated and have Ej > 25 GeV and |77 c i uster | < 2.47, exclud- 
ing the barrel-endcap transition region (1 .37 < |?7clusterl < 1.52), 
where Ej is the transverse energy and 77 c iuster is the pseudorapid- 
ity of the calorimeter energy cluster associated with the candi- 
date. The electron is defined as isolated if the Ej deposited in the 
calorimeter and not associated with the electron in a cone in j]-<p 
space of radius AR = 0.2 is less than 4 GeV. The muons must 
also be isolated and have pj > 20 GeV and \t]\ < 2.5. For iso- 
lated muons, both the corresponding Ej and the analogous track 
isolation transverse momentum (pj) must be less than 4 GeV 
in a cone of AR = 0.3. The track isolation pj is calculated 
from the sum of the track transverse momenta for tracks with 
pj > 1 GeV around the muon. Jets are reconstructed with the 
anti-A:, algorithm Il22ll with a radius parameter R = 0.4, starting 



"The fraction of events with m(*(- < 40 GeV is estimated to be less than 
0.2% of the total after all selections. 



from energy deposits (clusters) in the calorimeter reconstructed 
using the scale established for electromagnetic objects. These 
jets are then calibrated to the hadronic energy scale using pj- 
and ^-dependent correction factors obtained from simulation 
1 25]. The jet candidates are required to have pj > 25 GeV 
and \r]\ < 2.5. Jets identified as originating from a b quark (b- 
tag) by a vertex tagging algorithm are those that pass a decay 
length significance cut corresponding to an efficiency of 70% 
for fo-quark jets from it events and a 1% efficiency for light- 
quark and gluon jets J2, Hfl . 

The missing transverse momentum is constructed from the 
vector sum of all calorimeter cells with \rj\ < 4.5, projected 
onto the transverse plane. Its magnitude is denoted E 1 ™™. The 
hadronic energy scale is used for the energies of cells associated 
with jets; r candidates are treated as jets. Contributions from 
cells associated with electrons employ the electromagnetic en- 
ergy calibration. Contributions from the pj of muon tracks are 
included, removing the contributions of any calorimeter cells 
associated with the muon. 

4.1. t Reconstruction and Identification 

The reconstruction and identification of hadronically decay- 
ing t leptons proceeds as follows: 

1 . the r candidate reconstruction starts by considering each 
jet as a t candidate; 

2. energy clusters in the calorimeter associated with the r 
candidate are used to calculate kinematic quantities (such 
as Ej) and the associated tracks are found; 

3. identification variables are calculated from the tracking 
and calorimeter information; 

4. these variables are combined into multivariate discrimi- 
nants and the outputs of the discriminants are used to sep- 
arate jets and electrons misidentified as t leptons decaying 
hadronically from t leptons. 

Details are given in Ref. 9. In this analysis the outputs of BDT 
discriminants are used. 

Reconstructed r candidates are required to have 20 GeV < 
Ej < 100 GeV. They must also have \rj\ < 2.3, and one, two 
or three associated tracks. A track is associated with the t can- 
didate if it has pj > 1 GeV and is inside a cone of AR < 0.4 
around the jet axis. The associated track with highest pj must 
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have pi > 4 GeV. The charge is given by the sum of the charges 
of the associated tracks, and is required to be non-zero. The 
probability of misidentifying the r lepton charge sign is about 
1%. The charge misidentification rate for muons and electrons 
is negligible. 

If the t candidate overlaps with a muon (pj > 4 GeV, no iso- 
lation required) or an electron candidate within AR({, t) < 0.4, 
the t candidate is removed. To remove electrons misidentified 
as t leptons, an additional criterion is used that relies on a BDT 
trained to separate r leptons and electrons (BDT e ) using seven 
variables shown to be well modelled by comparing Z — > <? + e~ 
and Z — ♦ t + t~ events in data and in MC simulation. The vari- 
ables were chosen after ranking a large set by their effective- 
ness. 3 The most effective variables for BDT e are E/p, the 
EM fraction (the ratio of the t candidate energy measured in 
the EM calorimeter to the total r candidate energy measured in 
the calorimeter), and the cluster-based shower width. The BDT 
output tends to be near 1 (0) if the r candidate is a t lepton 
(electron). The r candidate is required to satisfy BDT e > 0.51; 
85% of reconstructed t leptons decaying hadronically satisfy 
that requirement as measured in Z — > t + t - events. The addi- 
tional rejection for electrons is a factor of 60. 

The majority of objects reconstructed as r candidates in a 
multi-jet environment are jets misidentified as r leptons (fake 
r). Another BDT (BDT,) based on eight variables is used 
to separate r leptons in t candidates with one track (denoted 
Ti) from such jets. For candidates with more than one track 
(denoted r?) BDT, includes ten variables. The most effective 
variables for BDT j are calorimeter and track isolation, cluster- 
based jet mass, and the fraction of energy within AR = 0. 1 of 
the jet axis. The BDT, distributions are fit with templates for 
background and signal to extract the number of r leptons in 
the sample. Details are given in Section [6] The fake r back- 
ground in the T3 sample is significantly higher than in the t\ 
sample, leading to very different BDT, distributions. Hence in- 
dependent measurements are carried out for T\ and T3 candidate 
events and the results are combined at the end. If there is a t\ 
and a T3 candidate in the event, the t\ candidate is kept as the 
probability that the t\ is a t lepton is much higher. If there are 
two T\ or T3 candidates, both are kept. 

4.2. Event Selection 

For this analysis, events are selected using a single-muon 
trigger with a pi threshold of 18 GeV or a single-electron trig- 
ger with a pj threshold of 20 GeV, rising to 22 GeV during 
periods of high instantaneous luminosity. The offline require- 
ments are based on data quality criteria and optimised using 
Monte Carlo simulation: 

• a primary vertex with at least five tracks, each with pj > 
400 MeV, associated with it; 



The effectiveness is quantified by quadratically summing over the change 
in the purity between the mother and daughter leaves for every node in which 
the given variable is used in a decision tree. 



• one and only one isolated high-px muon and no identified 
electrons for the p+r channel, or one and only one isolated 
electron and no isolated muons for the e + r channel; 

• at least one r candidate (as defined in Section 14. II ); 

• at least two jets not overlapping with a r candidate, i.e. 
Aft(r,jet) > 0.4; 

• E™ 1 ™ > 30 GeV to reduce the multi-jet background, and 
the scalar sum of the pj of the leptons (including r), jets, 
and E™ lss must be greater then 200 GeV, to reduce the 
W+jets background. 

The { + t samples are divided into events with no jets iden- 
tified as a fo-quark jet (0 b-tag control sample) and those with 
at least one such jet (> 1 fe-tag it sample). The fo-tag sample 
is used to estimate the background in the > 1 b-tag it sample. 
Each sample is split into two, one with the t candidate and I 
having the opposite sign charge (OS), and the other one with t 
and I having the same sign charge (SS). While the r candidates 
in the SS samples are almost all fake r leptons, the OS samples 
have a mixture of r leptons and fake t leptons. The numbers of 
observed and expected events in the above samples are shown 
in TableQ] All processes contribute more events to OS than SS 
because of the correlation between a leading-quark charge and 
the lepton charge, except for the multi-jet channel contribution 
which has equal number of OS and SS events within the uncer- 
tainties. The ^+jets entry includes the contribution from events 
with t leptons when the t candidate is actually a fake t. The 
t entries require the reconstructed t candidate be matched to a 
generated r lepton. The matching criterion is AR < 0. 1 between 
the t candidate and the observable component of the generated 
t lepton. 

To estimate the multi-jet background from data, an event se- 
lection identical to the p + r (e + t) event selection except for an 
inverted muon (electron) isolation cut is used to obtain a multi- 
jet template for the shape of the transverse mass, mj. 4 The 
normalization of each selected data sample is obtained by fitting 
the mj distribution of the selected data samples with the multi- 
jet template and the sum of non-multi-jet processes predicted 
by MC, allowing the amount of both to float. The uncertainty 
on the multi-jet background is estimated to be 30%. However, 
because of the subtraction method discussed in Section [5] the 
multi-jet background plays no role in the cross-section mea- 
surement. There are small differences between the total number 
of events predicted and observed which motivate using data as 
much as possible to estimate the background. 

As one can see from Table [1] the t leptons are almost all 
in the OS sample and come mainly from two sources: Z — > 
t + t ~, which is the dominant source in the sample with b- 
tag, and it — > I + t + X which is the dominant source in the 
sample > 1 &-tag. The sources of fake r leptons are also quite 
distinct between the b-tag and the > 1 b-tag samples: the 
first is mainly W/Z+jets with small contributions from other 
channels, the second is mainly it. 



4 m T = J(E£ + £™ ss ) 2 - (p f x + Ef' si ) 2 - (p r y + £™ ss ) 2 . 
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Table 1: Number of f+T candidates for Monte Carlo samples and data. ii({+e) are it events with one identified lepton and an electron reconstructed as a T candidate. 
ii(i + jets) are it events with one identified lepton and a jet reconstructed as a r candidate. ^+jets are events with one identified lepton and a jet reconstructed as a t 
candidates from sources other than it and multi-jets. Sources contributing to jet fakes are W+jets, Z+jets, single top-quark and diboson events. Wt(C + r) is W + t 
production with one W decaying to f and another to t. Excepting multi-jets the uncertainties are statistical only. MC samples are normalized to the data integrated 
luminosity 
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5. Background Models 

The jet origin can strongly influence the r-lepton fake prob- 
ability. Due to their narrow shower width and lower track mul- 
tiplicity, light-quark jets have a higher probability of faking a t 
lepton than other jet types. Thus the BDT ; distributions have 
a strong dependence on the jet type. It is therefore crucial to 
build a background model which properly reflects the jet com- 
position in order to correctly estimate the fake t contamination 
in the signal region. Deriving this background model from con- 
trol regions in data rather than MC simulation is preferable in 
order to avoid systematic effects related to jet composition in 
the MC models. 

The gluon component of the fake r leptons is charge symmet- 
ric; therefore it is expected to have the same shape in SS events 
as in OS events and should contribute the same number of fake 
t leptons in each sample. The contribution of fake t leptons 
from gluons can be removed by subtracting the distribution of 
any quantity for SS events from the corresponding distribution 
for OS events. The multi-jet background also cancels, as can 
be seen in Table [T] The resulting distributions are labeled OS- 
SS. Similarly, since each sample is expected to have an almost 
equal contribution from fo-jets and fe-jets, the small Z?-jet com- 
ponent should also be removed by OS-SS (asymmetric single b 
production is negligible compared to bb production). The only 
jet types remaining in the OS-SS distributions are light-quark 



jets. MC studies indicate that the BDT; distributions of c-quark 
jets misidentified as r leptons are not noticeably different from 
those of light-quark jets. 

One can construct a background BDT, distribution from the 
b-tag data sample by subtracting the expected amount of true 
r signal. The signal is mainly from Z — > t + t~ and can be 
reliably predicted from MC. A control sample dominated by 
W+jets events is considered as a check. The latter sample is 
selected by requiring events with a muon and a r candidate, no 
additionaljets, £™ ss > 30 GeVand40 GeV < m T < 100 GeV. 
According to MC simulation, in W+jets events where exactly 
one jet is required, 90% of the fake t leptons are from light- 
quark jets and 10% from gluons. This sample is labeled W + 1 
jet. 

The BDT; background shapes for the OS-SS b-tag and > 1 
b-tag data samples are not identical to the W + 1 jet distributions 
for two reasons: (1) the shape depends on the jet multiplicity, 
(2) different OS/SS ratios are observed in the samples. The 
dependence on the OS/SS ratio comes from the differences in 
jet fragmentation for a leading particle with the opposite charge 
and the same charge as the initial quark. MC studies of the ratio 
of OS-SS BDT, background distributions derived from W + 1 
jet and > 1 b-tag show that significant corrections are needed 
(30% for BDT, > 0.8, a region dominated by the true t signal). 
For the b-tag sample the corresponding corrections are much 
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smaller (5% in the same region). Both the b-tag and the W + 
1 jet data samples are used to obtain statistically independent 
estimates of the background in the > 1 b-tag sample. 

Two different approaches are used for constructing back- 
grounds in the > 1 b-tag data sample. One, used by the fit 
method (Section[6]l, is to reweight the BDTy distribution of the 
background bin-by-bin using the MC -based ratio of the > 1 b- 
tag background to the background model. In this case the 
b-tag sample is preferred as it requires smaller corrections de- 
rived from MC simulation; the W + 1 jet is used as a cross 
check. The other approach is to split the background into bins 
of some variable within which the shapes of BDT; distributions 
of the background model are close to those from the > 1 b-tag 
background. This approach, used in the Matrix Method cross 
check (Section fc.ll l. avoids using MC corrections, but assumes 
the data and MC simulation behave similarly as function of the 
binning variable. 



6. Fits to BDT y Distributions 

The contribution from t 1 — > i +t+X signal is derived from the 
> 1 b-tag data sample by ax 2 fit to the OS-SS BDT, distribu- 
tion with a background template and a signal template. The pa- 
rameters of the fit are the amount of background and the amount 
of signal. The shapes of the templates are fixed. 

Two background templates corrected by MC, as discussed in 
Section [5] are used: one derived from b-tag data, the other 
from the W + 1 jet data sample. The signal BDT, templates for 
b-tag and > 1 b-tag are derived from r leptons in it and Z — > 
t + t~ MC simulation. Contributions to the BDT; distributions 
from electrons passing the BDT e cut cannot be distinguished 
from t leptons so they are treated as part of the signal. 

The uncertainty on the background templates is determined 
by the numbers of data and MC simulated events. The signal 
template for the b-tag control sample also has non-negligible 
statistical uncertainty (2% for T\, 5% for t-j) because of the low 
acceptance. 

The fitting procedure was tested extensively with MC simu- 
lation before applying it to data. In the fits to the > 1 b-tag data, 
applying MC corrections to the b-tag background template in- 
creases the statistical uncertainty but raises the measured cross 
section by only 1%. 

Figure[T]shows the BDT; (OS-SS) distributions of (+t events 
with b-tag and the b-tag background template after subtract- 
ing the expected number of t leptons and applying the MC cor- 
rections. The t signal is mostly Z — > t + t~ events with a small 
contamination of electrons faking r leptons (from 1 1 — > I + e + X 
and Z — » e + e~) and a small contribution from it — > t + r + X. 
The uncertainty on the background template includes the sta- 
tistical uncertainty of the correction, the statistical uncertainty 
from MC and the b-tag data uncertainty. 

Figure |2] shows the result of the fit to the > 1 b-tag samples. 
The rlepton signal is mostly it — > C+t+X with a small contam- 
ination of misidentified electrons (estimated by applying fake 
probabilities derived from data), and small contributions from 
Z — > t + t~ events and single top-quark events (estimated from 
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Figure 1: BDT; (OS-SS) distributions of £+ t (e and ji combined) events in the 
fo-tag data (black points). The expected contributions from t and e are shown 
as a solid red line. The derived background templates are shown as dashed 
histogram with shaded/blue statistical uncertainty bands. The shapes of these 
background templates are used for the fits to the > 1 b-tag distributions after 
applying MC corrections. Top is for T\ , bottom for T3 . 



MC simulation). These contributions are subtracted from the 
number of signal events before calculating the cross section. 
The fit results using the background templates derived from 
b-tag data and W + 1 jet data are shown in Table [2] The results 
are consistent with each other within the statistical uncertainties 
of the background templates. The BDT, distributions for ti and 
T3 are fitted separately. The combined I + ti results are obtained 
by fitting the sum of the distributions. After adding ( + t\ and 
I + T3 signals obtained from a x 1 fit to the combined e + t and 
jj + T distributions and subtracting the small contributions to the 
signal from Z — > t + t~, Z — > e + e~ and tt — > e + ( (given in 
Table© the results are 840 ± 70 (243 ± 60) tt -> I + n (t 3 ) + X 
events. The uncertainty is from the fit only and does not in- 
clude systematic uncertainties. The results are in good agree- 
ment with the 780 + 50 (243 + 60) events obtained with the 
W + 1 jet background template and consistent with the number 
expected from MC simulation, 726 ± 19 (217 ± 10). Note that 
the fit uncertainty is dominated by the uncertainty on the back- 
ground template, thus the statistical uncertainties of the results 
with the two different background templates are not strongly 
correlated. 

Figure |3]shows the OS-SS distribution of the number of jets 
for > 1 b-tag events adding all channels for two BDTj regions: 
BDT, < 0.7, which is dominated by it -> I + jets, and BDT; > 
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Figure [4] shows the invariant mass of a selected jet with the t 
candidate for BDT, < 0.7 and BDT, > 0.7 for events with 
a r candidate and three or more jets. The selected jet is the 
highest pi untagged jet in events with more than one b-tag and 
the second highest pj untagged jet in events with one b-tag. 
The distribution shows clearly the presence of a W decaying 
to two jets in the BDT, < 0.7 region dominated by it — > t + 
jets. The mass distribution in the BDT, > 0.7 signal region 
is significantly broader as expected for it — > i + r + X. The 
signal and background shown in these figures are based on the 
fit using the b-tag background template. 
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Figure 2: BDTy (OS-SS) distributions of C + r in the > 1 i-tag sample. The 
normalisation of each template is derived from a fit to the data. The fitted 
contributions are shown as the light/red (signal), dashed/blue (background de- 
rived from i-tag after applying MC corrections) and dark/black (total) lines. 
Shaded/blue bands are the statistical uncertainty of the background template. 



Table 2: Results of template fits to fi + r, e + r and the combined BDT, distri- 
butions. The combined results are obtained by fitting the sum of the /i + t and 
e + t BDTj distributions. The first column gives the channel and the second 
the r type. The third column shows the extracted signal (sum of r leptons and 
electrons misidentified as r leptons) with the background template derived from 
fe-tag data distributions. The fourth column shows the extracted signal with 
the background template derived from W + 1 jet. The uncertainties are from 
the uncertainties in the fit parameters and do not include the systematic uncer- 
tainties. The MC columns give the expected r signal and the expected number 
of it — > / + r events after subtracting the contribution from non-?f events to the 
signal, assuming the theoretical it cross section (164 pb). 







Background template 


MC 








b-tag 


W + 1 jet 


Signal 


it 


p + T 


T\ 


490 ± 40 


456 ± 32 


432 


388 




T3 


135 ± 33 


130 + 50 


126 


116 


e + t 


Tl 


440 ±50 


430 ± 50 


388 


338 




T3 


116 ± 32 


120 ± 28 


114 


101 


Combined 


Tl 


930 ± 70 


860 ± 50 


820 


726 




1"3 


260 ± 60 


260 ± 40 


239 


217 



0.7, in which the largest contribution is from it — » ( + t + X. As 
expected, the multiplicity of jets peaks at four when BDTy < 
0.7 and three when BDT,- > 0.7 (the r is counted as a jet). 
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Figure 3: OS-SS number of jets distributions for events with at least one fo-tag. 
The fi + t and e + r channels have been summed together. The solid circles 
indicate data and the histograms indicate the expected signal and backgrounds. 
The normalisation of the expected signal and the backgrounds are based on the 
fit result. The uncertainty includes statistical and systematic contributions. The 
fraction of each background is estimated from MC. Top is for BDT ; < 0.7, 
bottom forBDTy > 0.7. 



6.1. Check with Matrix Method 

From Figures |3]and|4]one can see that a BDT ; > 0.7 require- 
ment separates well a region dominated by it — > f+jets from a 
region dominated by it — > I + t + X. One can extract the signal 
from the same OS-SS > 1 b-tag sample used by the fit method 
via a matrix method. All r candidates are labeled "loose", and 
t candidates with BDTj > 0.7 are labeled "tight". The proba- 
bility that the loose t candidates are also tight t candidates, for 
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Figure 4: OS-SS invariant mass of jet and r candidate for events with at least 
one b-tag. The jet is the highest pj untagged jet in events with more than one 
fe-tag and the second highest pj untagged jet in events with one />tag. The 
p + r and e + r channels have been summed together. The solid circles indicate 
data and the histograms indicate the expected signal and backgrounds. The 
normalisation of the expected signal and the backgrounds are based on the fit 
result. The uncertainty includes statistical and systematic contributions. The 
fraction of each background is estimated from MC. Top is for BDTy < 0.7, 
bottom for BDT; > 0.7. 
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both t leptons and fake r leptons, is defined as 



N 



tight 



N 



^real 



N 
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ffake 



tight 
fake 
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fake 



where the "real" subscript denotes r lepton, the "fake" subscript 
denotes fake r and N is the number of t candidates. The number 
of "tight" t leptons is then given by 



N 



tight 



real 



N 



tight 



£fake 



data 



^real ^fake 



data 



M"eal 



iV data >■ 



The value of efake is measured utilizing the OS-SS BDT, dis- 
tributions from the background control samples; e rea i is derived 
from MC and was tested using Z — > t + t~ events. This method 
uses the binning approach described in Section[5]to estimate the 
background. Values of ef a k e and e re ai are measured separately 
for three EM-fraction bins. The EM-fraction, the ratio of the 
energy measured in the EM calorimeter to the total t candidate 
energy measured in the calorimeter, is an effective variable for 
splitting the data into regions where the shapes of MC OS-SS 
BDTj distributions for the W+l jet background template and 



Table 3: Number of signal events obtained with the matrix method forjU+T, e+T 
and the combined channels. The first column gives the channel and the second 
the r type. The third column shows the extracted signal with the background 
template derived from b-tag data distributions. The fourth column shows the 
extracted signal with the background template derived from W+l jet. In order 
to compare the matrix method results to the fit results the number of signal 
events shown is 2 ^J^/e^ where e rea i is the e rea | averaged over the three 
EM-fraction bins. The uncertainties are statistical only. 







Background template 






b-tag 


W + 1 jet 


fX + T 


Tl 


460 ± 50 


440 ± 50 




T3 


130 + 40 


105 ± 35 


e + T 


Tl 


420 ± 60 


350 ± 50 




T3 


140 ± 40 


160 ± 40 


Combined 


Tl 


880 + 70 


800 ± 70 




T3 


270 ± 60 


260 ± 60 



the > 1 b-tag background are similar. Table [3] shows the num- 
ber of signal events obtained with the matrix method using the 
background derived from the b-tag data sample and from the 
W + 1 jet data sample. The numbers in each pair are in good 
agreement and consistent with the numbers obtained by fitting 
the OS-SS BDT ; distributions. 

6.2. Systematic Uncertainty 

Several experimental and theoretical sources of systematic 
uncertainty are considered. Lepton trigger, reconstruction and 
selection efficiencies are assessed by comparing the Z — > £ + i~ 
events selected with the same object criteria as used for the it 
analyses in data and MC. 

Scale factors are applied to MC samples when calculating 
acceptances to account for any differences between predicted 
and observed efficiencies. The scale factors are evaluated by 
comparing the observed efficiencies with those determined with 
simulated Z boson events. Systematic uncertainties on these 
scale factors are evaluated by varying the selection of events 
used in the efficiency measurements and by checking the stabil- 
ity of the measurements over the course of data taking. 

The modeling of the lepton momentum scale and resolution 
is studied using reconstructed invariant mass distributions of 
Z — > l + t~ candidate events and used to adjust the simulation 
accordingly II231 12411 . 

The jet energy scale (JES) and its uncertainty are derived 
by combining information from test-beam data, LHC collision 
data and simulation I25ll . For jets within the acceptance, the 
JES uncertainty varies in the range 4-8% as a function of jet 
pj and 77. Comparing MC and data the estimated systematic 
uncertainties are 10% and 1-2% for the jet energy resolution 
(JER) and the efficiency, respectively. The uncertainty on the 
efficiency of the ^-tagging algorithm has been estimated to be 
6% for Z?-quarkjets, based on Z?-tagging calibration studies [26]. 

The uncertainty in the kinematic distributions of the it sig- 
nal events gives rise to systematic uncertainties in the sig- 
nal acceptance, with contributions from the choice of genera- 
tor, the modeling of initial- and final-state radiation (ISR/FSR) 
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and the choice of the PDF set. The generator uncertainty is 
evaluated by comparing the MC@NLO predictions with those 
of POWHEG HI HI interfaced to either HERWIG or 
PYTHIA. The uncertainty due to ISR/FSR is evaluated using 
the AcerMC generator Q interfaced to the PYTHIA shower 
model, and by varying the parameters controlling ISR and FSR 
in a range consistent with experimental data 112 ill . Finally, the 
PDF uncertainty is evaluated using a range of current PDF 
sets H3 lL l32l 13311 . The dominant uncertainty in this category 
of systematic uncertainties is the modelling of ISR/FSR. 

The r ID uncertainty is derived from a template fit to a 
Z — > t + t~ data sample selected with the same p and r candi- 
date requirements as the sample for this analysis, but with fewer 
than two jets and mj < 20 GeV to remove W+jets events. The 
fit relies on the W + 1 jet data sample for a background template 
and Z — > r + r~ MC events for a signal template. The uncer- 
tainty includes the statistical uncertainty of the data samples, 
the uncertainty in the Z/y* cross section measured by ATLAS 
1 34] (excluding luminosity uncertainty) and jet energy scale un- 
certainty. It also includes the uncertainty on the number of 
misidentified electrons (< 0.5%, determined from Z — > e + e~ 
data). 



Table 4: Relative systematic uncertainties, in %, for the cross-section measure- 
ment. The first column gives the source of systematic uncertainty, ID/Trigger 
stands for the combined uncertainty of lepton identification and lepton trigger. 
The r ID uncertainty includes electrons misidentified as r leptons. The second 
and third columms give the channel. 



Table 5: Measured cross section from the T[ and T3 samples, as well as the com- 
bination (T1+T3) for each channel separately. The uncertainty in the integrated 
luminosity (3.7%) is not included. 



Source 


H + T 


e + t 


p (ID/Trigger) 


-1.1 /+1.5 




e (ID/Trigger) 




±2.9 


JES 


-2.0/+2.2 


-1.9 /+2.8 


JER 


+ 1.0 


±1.2 


ISR/FSR 


+4.8 


±3.5 


Generator 


+0.7 


±0.7 


PDF 


+2.0 


±2.1 


fo-tag 


-7.7/+9.0 


-7.5/+8.9 


T] ID 


-3.0/+3.2 


-2.7/+3.0 


T 3 ID 


-3.1/+3.4 


-2.9/+3.2 



The effect of these variations on the final result is evaluated 
by varying each source of systematic uncertainty by ±lcr, ap- 
plying the selection cuts and recalculating the cross section. 

The uncertainties obtained for the fit method using the b- 
tag background template are shown in Table |4] The systematic 
uncertainties for the matrix method are very similar. The uncer- 
tainty on the measured integrated luminosity is 3.7% 113311 . 



7. Measuring the tt Cross Section 

The cross section is derived from the number of observed 
OS-SS signal events in the > 1 fo-tag data sample assuming the 
only top quark decay mode is t — > Wb, and subtracting from 
that number the small contribution from it — > e + 1 (from elec- 
trons faking r leptons) and t leptons from Z — > t + t~ (TableQ}. 





ii A- T 

LA 1 1 


Tl 


189 ± 16 (stat.) ± 20 (syst.) pb 


T3 


180 ±40 (stat.) ±21 (syst.)pb 


TI+T3 


186 ± 15 (stat.) ± 20 (syst.) pb 




e + t 


Tl 


190 ± 20 (stat.) ± 20 (syst.) pb 


T3 


170 ±50 (stat.) ±21 (syst.)pb 


T1+T3 


187 ± 18 (stat.) ±20 (syst.)pb 



The systematic uncertainties are estimated as the quadratic sum 
of all uncertainties given in Table |4] which includes the uncer- 
tainty from the subtraction. 

The results are given separately for ti and T3 and then com- 
bined (weighted by their statistical uncertainty and assuming 
all systematic uncertainties other than from r ID are fully cor- 
related). The results using the b-tag background template are 
shown in Table [5] 

The results for the p + T and e + r channels are combined tak- 
ing into account the correlated uncertainties using the BLUE 
(Best Linear Unbiased Estimator) technique 03611 . Combining 
them does not improve the systematic uncertainty as the sys- 
tematic uncertainties are almost 100% correlated. 

The results for each lepton type are: 

p + r: cr fl = 186 ± 15 (stat.) ± 20 (syst.) ± 7 (lumi.) pb, 

e + r: cr, 7 = 187 ± 18 (stat.) ± 20 (syst.) ± 7 (lumi.) pb, 
Combining both channels one obtains: 

cr (7 = 186 ± 13 (stat.) ± 20 (syst.) ± 7 (lumi.) pb 

To check the fit measurements, the cross sections can be cal- 
culated using the matrix method and the results obtained with 
the W + 1 jet background to minimize the correlation with the 
fit results. The combination of the matrix method and the fit 
results with the BLUE method show they are compatible at the 
45% and 10% confidence level for p + r and e + r, respectively. 

8. Conclusions 

The cross section for ti production in pp collisions at 7 TeV 
has been measured in the yU + r and the e + r channels in which 
the r decays hadronically. The number of r leptons in these 
channels has been extracted using multivariate discriminators 
to separate r leptons from electrons and jets misidentified as 
hadronically decaying r leptons. These numbers were obtained 
by fitting the discriminator outputs and checked with a matrix 
method. Combining the measurements from p + t and e + r 
events, the cross section is measured to be 

cr, 7 = 186 ± 13 (stat.) ± 20 (syst.) ± 7 (lumi.) pb, 
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in good agreement with the cross section measured by ATLAS 
in other channels with the cross-section measurement 

by the CMS Collaboration JUlltl and with the SM prediction, 
164+1' pb H. 
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